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Abstract 

We extend the recent D = 5 results of Dias, Horowitz and Santos by finding asymptotically AdS 
rotating black hole and boson star solutions with scalar hair in arbitrary odd spacetime dimension. 
Both the black holes and the boson stars are invariant under a single Killing vector field which co- 
rotates with the scalar field and, in the black hole case, is tangent to the generator of the horizon. 
Furthermore, we explicitly construct boson star and small black hole (r-|_ ^ £) solutions perturbatively 
assuming a small amplitude for the scalar field, resulting in solutions valid for low energies and angular 
momenta. We find that just as in I? = 5, the angular momentum is primarily carried by the scalar 
field in D > 5, whereas unlike D = b the energy is also primarily carried by the scalar field in D > 5; 
the thermodynamics in D = 5 are governed by both the black hole and scalar field whereas in D > 5 
they are governed primarily by the scalar field alone. We focus on cataloguing these solutions for the 
spacetime dimensions of interest in string theory, namely _D = 5, 7, 9, 11. 
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1 Introduction 



In general, finding analytic solutions to Einstein's equations is not an easy task; a common strategy is to 
assume a set of symmetries and input a suitable ansatz into the field equations. This is not guaranteed 
a priori to produce a consistent solution but there exist theorems on spacetime structure which aid in 
choosing an appropriate ansatz. For instance, it is well-understood that in any number of dimensions a 
stationary spacetime must also be axisymmetric [21 El S] , that four-dimensional black holes can only have 
a spherical topology [S], or that black holes in spacetimes with non-negative cosmological constant have 
no hair [6l [7] when coupled to "ordinary" Maxwell- type matter fields. 

Static black holes without rotational symmetry were constructed perturbatively in situations where 
the horizon radius is close to the critical radius for instability of the Reissner-Nordstrom solution [1]. 
Essentially a magnetically charged Reissner-Nordstrom black hole whose horizon radius is less than a 
critical value near its inverse mass is unstable classically against the development of a nonzero vector 
meson field just outside the horizon. Such black holes have vector meson hair. For spacetimes with a 
negative cosmological constant the no-hair theorem no longer applies because the reflecting boundary 
conditions of AdS can support a non-trivial matter field. Despite this, all known asymptotically AdS 
black hole solutions had at least two Killing vectors, regardless of the presence of hair. Recently a novel 
solution was constructed [8] that considered five-dimensional Einstein gravity with negative cosmological 
constant minimally coupled to a massless complex doublet scalar field, describing "lumpy" scalar hair co- 
rotating with a black hole. This configuration is motivated by superradiance: scalar fields in AdS spaces 
can increase their amplitude by scattering off the horizon of a rotating black hole, which is then reflected 
back to the horizon by the AdS boundary conditions, leading to a further increase in amplitude. This 
process mines rotational energy from the black hole and the end result is lumpy scalar hair co-rotating at 
the same angular velocity with the black hole. The spacetime and scalar fields are collectively invariant 
under a single Killing vector, which is tangent to the generator of the horizon. 

In this paper, we extend the results of [8] to arbitrary odd spacetime dimension D > 5 and construct 
analytic boson star and black hole solutions perturbatively in the dimensionless scalar amplitude parame- 
ter e ^ 1. We do so by a cohomogeneity-1 ansatz, i.e. metric functions of the radial coordinate only, and 
a judicious choice of scalar fields whose stress tensor shares the symmetries of the metric. This ensures 
that the resulting equations of motion form a set of coupled ODEs instead of a system of coupled PDEs. 
We catalogue the results for the spacetime dimensions of interest in string theory, D = 5,7,9,11. We 
include the D = 5 results in the interest of having our paper be self-contained and because we find various 
discrepancies between our results and the results of [8] . These discrepancies are explicitly discussed when 
they appear in ^4.2.21 and ^4.2.3[ 

The rest of the paper is outlined as follows. In ^ we introduce the metric and scalar field ansatz 
and give the set of ODEs following from the equations of motion. Next, in ^ we construct perturbative 
solutions for the boson star in AdS using the scalar field amplitude as an expansion parameter. Pertur- 
bative small black hole solutions follow in ^by introducing a second expansion parameter, r^/£ <C 1. 
It is interesting to note that this two-parameter class of black hole solutions is related to both the boson 
star {r^ — t- 0) and Myers- Perry (e — t- 0) AdS solutions [9]. The thermodynamic properties of both the 
boson star and black hole solutions are discussed in ^ Finally, we conclude in ^with a discussion of 
our results and of future work needed in this area. 

2 Setup 

In this section we introduce the model for constructing hairy black holes and boson stars in arbitrary 
odd dimension and give the resulting equations of motion. 
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2.1 Metric and Scalar Field Ansatz 



We begin with D = n + 2 dimensional Einstein gravity with negative cosmological constant minimahy 
coupled to an ^^^y^-tuplet complex scalar field 
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(2.1) 



where we take the usual convention A = — "^^^'a . In order to obtain the desired symmetries in our 
solution, namely that the matter stress tensor has the same symmetries as the metric, we will need to 
Hopf fibrate our n-sphere. Thus, we consider only odd dimensions with n > 3 and propose the metric 
and scalar field ansatz 



-f{r)g{r)dt^ + 
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where Zi are complex coordinates such that dzidzi is the metric of a unit n— sphere. An explicit and 
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convenient choice for the Zi is 
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in which case dzidzi = {dx + Aidx^)'^ + gijdx^dx^ is the Hopf fibration of the unit n— sphere where 
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(2.5) 



and Qij is the metric on a unit CP 2 . In these coordinates the scalar fields are manifestly single- valued 
on the spacetime since x ^^'i 4'i have period 27r while the 6i have period ^. To verify our ansatz at this 



we 



point, we note that if we choose n = 3 and perform the coordinate transformation X = "0 ~ f ' ^ ~ 
recover exactly the ansatz considered in Ref. [8j. 

The form of the scalar fields is crucial to this construction and was first considered in jl5j : it is 

n-\-l 

clear from Eq. ()2.3p that the scalar fields can be viewed as coordinates on C^~. Given that n(r) is a 
function of r only, for each value of r, 11 traces out a round n-sphere with a time- varying but otherwise 
constant phase. On the other hand, constant r surfaces in the metric ()2.2p correspond to squashed 
rotating n-spheres. The stress tensor for the scalar field takes the form 



ab 



dJi*di,Yl + dJldbU* ]-gab[ djld'^li 



(2.6) 



which has the same symmetries as the metric (|2.2p since the first term is the pull-back of the round 
metric of the n-sphere and the second term is proportional to Qah- 
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Although the matter stress tensor has the same symmetries as the metric, the scalar fields themselves 
do not. Indeed, the metric ()2.2p is invariant under dt, as well as the rotations of CP~2~ while the 
scalar field (|2.3|) is only invariant under the combination 

K = dt + ujd^. (2.7) 

Therefore, any solution with non-trivial scalar field will only be invariant under the single Killing vector 
field given by (pT]) . 



2.2 Equations of Motion 

The equations of motion resulting from the action (|2.ip are Gab — "^2^V'* gafe = Tab and V^H = which 
ought to have non-trivial solutions by virtue of the matter stress tensor possessing the same symmetries 
as the metric. Indeed, inserting the ansatz (j2.2p and (j2.3p into the equations of motion yields a system 
of five coupled second order ODEs 



f - ^ [ ^ - T + -] + + T- + + + 



fr \ 6 6 J r i'^ r fg 

t2 /i , (n-l), 



2 4n^^l + i^/.j 2(n-3) 



(2.8) 



4H g' 1\ ^ {Er'^Vhy i^h^ -U^ 3(n + l) (n - 3)E 



(2.9) 
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,,, h' 2h' f frh'\ h^Pfl'^ 4(l-/i) / 2 (n + 1) \ , , 

^„ 4n2 ff 2h' 2H 2n + l\ , , 

n" _ ^nV^ _ / V n(^ - 17)^ _ {i + {n-i)h)u _ 

where H = /i + 11^ — ^^jl — ^'"^^V'' and a ' denotes differentiation with respect to r. In addition to these 
second order ODEs, the Einstein equations further impose two first order ODEs in the form of constraint 
equations, Ci = and C2 = 0. Explicitly, these are 



(/25^2{n-l)V 



C.= ''Z.n^s ' +^9hE (2.13) 



, PhQ'^ E{r-+^hy (hfYh' n{fhr-~'y -U^ n+1 

^ pg 4fg fhr^+^ 4fh^ 2fhr^ fhr^ 2fP ' ^ ' ' 

We note that inserting n = 3 into the above equations of motion yields the 5D results of [8]. Addi- 
tionally, it is interesting to note the presence of terms proportional to n — 3 in Eqs. ()2.8p and ()2.9p . which 
are therefore absent in 5 dimensions; one might expect such terms to change the physics of the higher 
dimensional solutions, though at the perturbative level this may not be apparent. Before we move on, 
we emphasize that the above equations of motion are exact for arbitrary odd dimension D > 5; we have 
explicitly verified this for D = 5, 7, 9, 11. 
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3 Perturbative Boson Stars 



In this section, we present the boundary conditions that define a boson star and then use these to 
construct such solutions as perturbations around AdS. The boson star is a horizonless solution for the 
matter configuration we are using and, furthermore, can be viewed as a warm-up problem to the hairy 
black hole solutions we will later construct. The expansion is carried out in orders of the scalar field 
condensate parameter e and we give results up to order for the spacetime dimensions of interest in 
string theory, namely D = 5,7, 9, 11. As a perturbative construction, these results will only be valid for 
small energies and angular momenta. 



3.1 Boson Star Boundary Conditions 
Boson Star Origin 

Boson stars are smooth, horizonless geometries, which means that all metric functions must be regular 
at the origin. Furthermore, due to the slow physical rotation of points as r — ?• 0, surfaces of constant 
t in the vicinity of the origin ought to be described by round n-spheres with r being the proper radial 
distance. To find the boundary condition on IT, we multiply (|2.12p by and note that 11 must vanish at 
the origin in order to yield consistent equations of motion. Thus, the boundary conditions at the boson 
star origin take the form 

fUo = ^ + 0{r^), gUo = 9{0) + 0{r), hU, = l + 0{r% (3.1) 



Boson Star Asymptotics 

In order to simplify the asymptotic boundary conditions, we first make note of a residual gauge freedom. 
It is straightforward to show that the transformation 

X^X + >^t, 0-;>$7 + A, to ^ uj - X (3.2) 

for some arbitrary constant A, leaves both the metric (|2.2p and scalar field (|2.3p unchanged. We will use 
this gauge invariance to pick a frame which is not rotating at infinity, i.e. we use it to set — t- in the 
limit r — 7- oo. 

In the r — )• oo limit the boundary conditions for the boson star will asymptote to AdS with corrections 
for mass and angular momentum. To ensure the solution has a Newtonian potential of the correct strength 
for a spherically symmetric mass distribution, we impose an r"^""^) correction for /. Next, requiring 
our solutions to have finite masses and angular momenta means we must impose an r"^""'"^) fall-off for 
the corrections to g,h and 0. These considerations determine the boundary conditons for f,g,h, and 
up to some constants Cf,Ch, and Cq. The remaining boundary condition is set by requiring 11 to be 
normalizable, which means it must decay like r~^'^'^^\ Explicitly, the asymptotic boundary conditions 
are given by 

2 pn—1 (-1 pn+1 

hU^ = 1 + ^ + = ^ + 0(.-("+^)), (3.3) 

Here and in what follows, e provides a dimensionless measure of the amplitude of the scalar field. 

With the appropriate boundary conditions and the equations of motion at hand, we are now ready 
to construct perturbative boson stars in arbitrary odd dimensions. 
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3.2 Perturbative Boson Star 



We start by expanding our fields in terms of the scalar field amplitude, e, as follows: 

mm m 

F{r,e) =Y^F2,{r)e^' U{r,e) = Y,f^2^+l{r)e^'^' ^(e) = Y.Ci;2^e''' (3.4) 

i=0 1=0 i=0 

where F = {f,g,h,Q} is shorthand for each of the metric functions in (|2.2p . The metric functions 
are expanded in even powers of e while the scalar fields are expanded in odd powers. This allows a 
perturbative expansion as follows: start with global AdS at m = 0, then introduce non-trivial scalar 
fields into the AdS background, without back-reacting on the metric, by solving ()2.12p . The full set of 
equations of motion will then be satisfied up to order e. At next order, m = 1, the scalar fields due to 
ni(r) then source corrections to the gravitational fields F2{r), and these in turn back-react on the scalar 
fields n3(r). The equations of motion will then be satisfied up to order e^. The perturbative solution 
can, in principle, be obtained by this bootstrapping procedure up to arbitrary order, m. Note that we 
also must expand the frequency in even powers of e. This is because at the linear order, the frequency is 
determined by the scalar field alone but at the next order, the scalar field then back reacts on the metric 
inducing non-trivial frame-dragging effects which in turn have a non-trivial effect on the rotation of the 
scalar field. In practice, these corrections to co are found by imposing the boundary conditions. 
Global AdS is given by 

fo = l + -, 50 = 1, ho = l, no = (3.5) 

In this background, the most general massless scalar field solution to ()2.12p which is consistent with the 
asymptotic boundary conditions (|3.3p is given by 



ni(r) = ^2^1 



n + 2-uj£ n + 2 + uji n + 3 



J.2 _j_ 



(3.6) 



where 2F1 is the hypergeometric function. Now in order to satisfy the boundary conditions at the origin 
()3.ip we must further restrict u to 

u!i = n + 2 + 2k, /c = 0,1,2,... (3.7) 

where the non-negative integer, k, describes the various possible radial modes of the scalar field. Although 
in principle, any radial profile can be built up out of a linear combination of the radial modes, this 
introduces multiple frequency parameters, w^. This is inconsistent with the existence of the Killing 
vector field (j2.7p . Furthermore, the mode A; = yields the ground state while higher modes represent 
excited states. In what follows, we therefore choose = as the only mode present, in which case p.6p 
simplifies to 

ni(0 = -— --ii^- (3.8) 

(r2 + 2 

Proceeding up the perturbative ladder, we insert (|3.8p and the expansion ()3.4p into the equations of 
motion, expand in e and solve for order e^. In general, the solutions contain two constants of integration, 
which are then uniquely fixed by the boundary conditions. These fields, F2{r) are then inserted into the 
equation of motion for n(r) to find the correction to the scalar fields. This process can, in principle, 
be taken to arbitrary order in e. In practice, however, the expressions become rather unwieldy at higher 
orders making this increasingly difficult to accomplish. Up to order e^, we find the general solutions to 
take the form 



fir) = 1 + ^ - ^TT^^' - AS:, ^ + 0{e^) (3.9) 
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Hr) = l + p^^^,e^ + Oie^) (3.11) 

(r2+£2)— (^2_^£2)^- _^ £2)^3-^ 

where the fields {fn;s,o, gn\s,o, hn;s,o,^n;s,o,^n;s,o} are simple polynomials in r; in this notation s labels 
the order in e and n = D — 2 labels the spacetime dimensioE0- These fields are catalogued in Appendix B 
for n = 3, 5, 7, 9 up to order e^. Similarly, the explicit corrections to ujn for the boson star can be found 
by taking the r+ — )• limit of the general expressions for cj^ in Appendix A. 



4 Perturbative Black Holes 



We now turn our attention to the more interesting hairy black hole solutions. The presence of a horizon at 
r = r+ provides us with an additional length scale, or rather, another perturbative expansion parameter, 
r_|_/£, on top of the scalar field amplitude, e, making this a two-parameter family of solutions. To 
find the black hole solution then, we will perform a double expansion in the scalar field condensate 
parameter e and dimensionless horizon radius r+ / i. As pointed out in [8] , inserting this double expansion 
into the equations of motion results in a set of differential equations that cannot be simultaneously 
solved everywhere in the spacetime. This problem is circumvented by applying a matched asymptotic 
expansion; this procedure will be made explicit in ^4.2.31 Basically, this involves splitting the spacetime 
into two regions, a far-region (r S> r+) and a near-region (r_|_ < r <^ i), and solving the equations 
of motion in each separately. At each order in our solution we will have two arbitrary constants of 
integration. The boundary conditions can be used to fix one of these for each region, i.e. we can apply 
the asymptotic boundary condition in the far-region and the horizon boundary condition in the near- 
region. The remaining constant can be fixed by matching the solutions where the two regions overlap: 

<^ r <^ i. This admits a unique solution valid in the entire spacetime that satisfies both boundary 
conditions. This analysis is carried out in the following section and we give results up to C'(e'"(r_|_/£)") 
where m + n < 6 for the spacetime dimensions D = 5,7, 9, 11. As in the boson star case we will consider 
only the ground state hairy black holes with regard to the excitation of the scalar field and, again, these 
results will only be valid for small energies and angular momenta. 

The double expansion of our fields can be interpreted as placing a black hole inside a rotating boson 
star, or alternatively as placing a non-trivial scalar field around a small rotating black hole. Consequently, 
in the limit r-|_ — )■ we should recover the boson star of the previous section and similarly, in the limit 
e — )• we should recover an asymptotically AdS rotating black hole, i.e. a Myers-Perry-AdS black hole 
[SI [Tot \TT\ I12j . We then see that the two-parameter family of hairy black hole solutions is connected 
to the two-parameter family of Myers-Perry-AdS solutions. For the hairy black hole, the frequency is 
uniquely determined in terms of r_|_ and e, so in the e — )■ limit, the hairy black hole joins with a 
one-parameter subset of Myers-Perry-AdS black holes, whose horizon angular velocity is identified with 
the e — 7- limit of a;. To phrase this in terms of the space of solutions, the one-parameter family of 
boson stars corresponds to a line in the (e, a;)-plane, the Myers-Perry-AdS black holes correspond to the 
{u},r^)-p\ane, while the hairy black holes correspond to a sheet through the {e,u},r^) octant. This sheet 
meets the {e,u}) plane on the line defining boson stars and it intersects the (w,r_|_) plane, which is where 

^The purpose of the "0" index on these coefficients will be clear when we compare to the perturbative black hole solutions 
which involves a second pertubation parameter. 
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it joins with the Myers-Perry- AdS family of solutions. This is represented schematically in the following 
Table. 



Order 




e2 g4 




(^)° 


Global AdS 


— )• Perturbative Boson Star — )• 


Exact Boson Star 




; 

Pert. MP- AdS BH 
i 


\ Perturbative \ 
Hairy 
\ Black Hole \ 










Exact MP- AdS BH 




Exact Hairy BH 



Before continuing, we explicitly present the Myers-Perry- AdS solution with equal angular momenta 
in the two planes of rotation. In terms of the metric ansatz ()2.2p . the metric functions take the form [8] 



/=7T + 1 



n-l 



1-^ + 



h = l + 



1 



where vm and a are related to the outer horizon, r+, and the angular velocity of the horizon, Qh = f^, as 



n = (4.1) 



M 



+ ^2- 



(4.2) 



Small black holes are then obtained by Taylor expanding the metric functions ()4.ip in powers of r_|_ /i. 



4.1 Black Hole Boundary Conditions 
Black Hole Horizon 

We wish to find non-extremal black hole solutions with scalar hair, so we begin by defining the location 
of the non-degenerate outermost horizon to be r = r_|.. Consequently, / must have a simple zero at r+ 
while all other metric functions must remain finite and non-zero. For the scalar field, one might expect 
H to vanish in the vicinity of the horizon since it was shown in |13[ [H] that one cannot have black holes 
inside boson stars. However this applies to static black holes and it was pointed out in [8] that this 
prohibition is removed if the black hole and the boson star are co-rotating. Indeed, if we require n(r+) 
to be finite and non-zero, then multiplying Eq. ()2.12p by shows that the equations of motion remain 
consistent across the horizon provided ^^(^4.) = oj. The boundary conditions at the black hole horizon 
are thus 

flr^r^ = 0{r - r+), g\ = g{r+) + 0{r - r+), h\ = h{r+) + 0{r - r+), 

, , (4.3) 

n\ =uj + 0(r-r+), H =0(r-r+). 



Black Hole Asymptotics 

The asymptotic boundary conditions for the black hole will be identical to those of the boson star since 
both are globally AdS with next-to-leading order terms accounting for mass and angular momentum. As 
in the case of the boson star, we are also free to exploit the gauge freedom in Eq. (j3.2p and choose to work 
in a frame which is non-rotating at infinity. Thus we will apply to the black hole the same asymptotic 
boundary conditions that we saw for the boson star in Eq. ()3.3p . 
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4.2 Perturbative Hairy Black Hole 



4.2.1 Far Region 

As discussed above, we start by performing a double expansion of our fields in e and r^ji as follows: 



2i 



n 



out 



(r,e,r+) = ^^n: 



out 

2j+l,2jV 



( 1^ ) (4.4) 



2i 



i=0 j=0 



i=0 j=0 



where F*"** 
As discussed in 



{/""*, (7°^*, 0°"*} is shorthand for each of the metric functions in the far-region ()2.2 
the scalar frequency a; has a similar expansion 



i=0 j=0 



2i 



(4.5) 



where this expansion holds in both the outer and inner regions since w is a globally defined constant. 

The perturbative expansion then proceeds similar to the boson star case except that instead of our 
background being global AdS we take it to be the Myers-Perry AdS black hole (|4.ip expanded in powers 
of r_(_/^ ^ 1. We can therefore immediately write down the field coefficients -FqIK'") t° arbitrary order 
while trivially satisfying the field equations and asymptotic boundary condition. With the complete r^^jl 
expansion in hand, the goal then is to introduce a non-trivial scalar field and solve the field equations 
in powers of e. As before, this is done by first inserting the fields F^'^^ir^ into Eq. ()2.12p and solving 

the order e equation to determine n^'2i order by order in r_|_/£. The presence of the scalar field then 
sources the gravitational fields at order e^, which are then determined by solving Eqs. (|2.8|) - (|2.11|) for 
^2 2i('^) order by order in r+/£. At the next order in e, there is a back-reaction on the scalar field, which 
is then similarly solved for from Eq. (|2.12|) order by order in r^ji. This iteration process continues with 
gravitational field corrections at every even order in e and scalar field corrections at every odd order in 
e. Due to the powers of that appear in the Myers-Perry solution, we carry out this procedure up to 
C(e°(r+/^)"'+^), C'(e^(r+/£)"+-^) and C'(e^(r+/£)""^). Applying the asymptotic boundary conditions we 
find the fields to have the structure: 



out I 



^2 

£2+1 r 



((n + 2)2 + l)r-ir!^+i 



n— 1 lln—\ 



r,n—l 



+ 



+ 



pn—1 fout 

^ /n;2,2 ^+ 

^n-3(^2 _j_ £2)71+2 £n-l 



+ 




+ e 



£n+l ' ^ I £n+3 I 



n;2,0 



(r2 + ^2) 



2n+3 



l + O 



r.n+3 



£n+3 



n(r2 + £2)n+2 



+ 



out 
i/n;2,2 



^n-3(^2 _|_ £2)n+3 £n-l 




(4.6) 



+ 



(4.7) 



/i°"*(r) 



(7.2 _|_ £2)2n+4 



l + O 



£n+3 

(n + 2)r+^ rl+^ 

^n+l fn+l 




(^2 _|_ £2)2n+3 



+ 




+ 



n;2,2 



„n-l 



(7-2 _|_ £2)n+l ' yn-3(j.2 _|_ £2)n+2 £n- 



+ 




(j2 _j_ £2)2n+3 



+ 
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n°"*(r) 



re 



m+l 



+ 



pn+ljiout 

^ ^^n;l,n-l 



„n-l 



J ^^n.;l,n+l ^+ 



r-r+3n„.3n 



+ 



TTOUt 

-^-'■n;3,n-l 



„n-l 



(r^+£^) 2 ^n-2(-^2 _^ ^2) — 2 — <^ 



+ 






(4.10) 



5(n+l) 

_(r2+£2)^^ 



are new fields which 



where the fields -Fn;s,o are the boson star fields of the previous section and 
enter at non-zero order in r+/£. In solving for these fields, one of the arbitrary constants from each 
second order ODE is fixed by the asymptotic boundary conditions while the other constant is fixed by 
the matching condition. The double expansion of w, which is also obtained by matching the inner and 
outer solutions, is catalogued in Appendix A and the fields with the matching condition already imposed 
are catalogued in Appendix B. We postpone a discussion of the matching procedure to ^4.2.31 Note that 
taking the r+ — t- limit of the above yields the boson star fields of ^3.21 as it ought to. 



4.2.2 Near Region 

Since we are constructing solutions perturbatively we are assuming low energies and angular momenta, 
so by construction r_|. <^ £. Just as in our asymptotic expansion we expanded in the dimensionless 
parameter e ^ 1, near the horizon we expand in the dimensionless parameter r^/i <^ 1. To accomplish 
this, we switch to the radial coordinate z = ir/r^ such that the horizon is located at z = i. This ensures 
that z/i > 1 is always large with respect to r^/i <^ 1 so we can safely expand our fields in powers of r^/i 
without needing to worry about competing effects at the same order. Note that this was also achieved 
in our asymptotic expansion since r/£ ^ 1 and the expansion parameter was e ^ 1. Thus, we perform 
our double expansion in the inner region using z as our radial coordinate. 
The gravitational and scalar field expansions are given by 



i=0 j=0 



2« ~. „. , /ri \ 2i+l 



n"'(.',e,r+) 



1=0 j=0 



2j+l,2i+l 



iz)e 



2j+i r_ + 



(4.11) 

(7™, /i*", fi*"} are the metric functions in the near-region. Recall, however, that the 



where F*" = {/ 
frequency uj is still given by Eq. (|4.5p since it is defined globally. 

Again, we start at order e*^ with the Myers- Perry- AdS solution, with gravitational fields FQ2i- We 
insert this into the equations of motion, appropriately transformed to equations of z, add a non-trivial 
scalar field by expanding to order e and solve the equations order by order in r^/i. Once all orders 
in r-^/i have been calculated up to the desired cutoff, the matching conditions must be imposed before 
continuing up the perturbative ladder. A discussion of this procedure is postponed to ^4.2.31 Next, the 
scalar field at order e sources the gravitational fields at order e^. These have to be solved order by order 
in r_|_/^ and matched to the far region expansions before the process is continued. The result of this 
procedure, with the matching conditions already imposed, up to 0{e"'{rj^/£)^) such that a + 6 < 6, is 
given by 



yU-l ~^ 1 £2 



pn-1 



pri+l ^ ^2 



+ ( 2(n + 2)w„;o,2^ + (n + 2f {{n + 2f - l] 



+1 



(4.12) 
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pn+l 

l-(n + 2)2' 



+ <^ - 2(n + 2)a;„;o,2^ + (n + 2)^ ((n + 2f - l) 



+ (n + 2) 



2 , rl /4r4fe] 



+ 



X ( [2x 



(4.13) 
1]) ) dx 



£4 



£3n+7 ^^1^2 



1 + (- + + (2(- + 2)-n;oy + (n + 2)2((n + 2)^ " l) j + O j 



2 / 4 

^n;2,2^ + 0^ ^ 



£4 



(4.14) 



(- + 2)^ + |^-„;0,2^ + (n + 2)3(z"+i - r+i)^| ^ + |^a;„,;o,4^ 



+ 



n+l 



r+^) (n + 2)2 ^ (^3£.„;o,2 + (n + 2) ((n + 2)2 - 1) j + (n + 2)^^-3;^ j | ^ 



(4.15) 
4 



^0Jn;2fi + (^^^n;2,o(0) - ^Wn;2,o) ( 1 ^TTT 1 + -^\^n-a,2 + 



-'^n;2,2 
rpn+2 



4(n + 2) 



yn+2^yn~l _ 
(^2_i) n+l 

4(yn-l _ 1) ■lri,2,2\ ) 



^n-X2^V 



2n 



n + 3 



,,"•+1 



y r2 



g||(P^[2,--l]) 



«n+l\ /^2 
+ (^f^n;4,0(0) - £^.;4,0) ( 1 " ^ ) + ^ ( ^ 



tf"(z) = e 



P 



22^^-1 



-"■n;l,3-'^- 



2£^ 

£2 



+ - ^^^iWX 



ri-l 



(--i)r[i^ 

P_^ r22;"-i - ll sJx) + nP_!L_ \2x'^-'^ - 1 
/ P_^ r2x"-i - ll : 

/-I n-1 / n-3 \ ^ 



(4.16) 

Snix) 



-dx 



n-1 



■2z"-i 



- 1 



rz/i 

J [2^"-' - 1 



P_^ r2x"~i - ll sJx) + nP^ r2x"-i - ll sJx 

71-1 '- -' '- -' 



■dx 



+ 6^ 



n — 1 



n-1 



- 1 



where the constants 2 ^^^"^ -^rri 3 calculated in ^4.2.31 and the constants -^^^3 1 hi Eq. 
()4.16p are given by 

r^u _55vr n _ WGJT^ _ 6403r2[£ 

^3;3,l - ' ^5;3,1 " ggp^ , ^7;3,1 " 14n2r[|] ' 

for n = 3, 5, 7, 9 respectively. The functions Pu[y] and Qu[y] are the Legendre functions of the first and 



n 110460ir2 I 

= 131 4.17 

^'^'^ 3175200r f ^ ^ 
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second kind, respectively, the functions s„(a;) are given in Appendix B and 



f n—3 \ / n—3 

2 \ / 2 



sn{^) = U + l)x"-3-2^J + (n + 2)2 |^ (j + l)^^^^ . (4.18) 

There are a few discrepancies with Ref. [8] (the n = 3 case) that warrant mentioning. In the 
order e(r+/£)^ term of Eq. ()4.16p . we used the matching condition outhned in the next subsection and 
determined that the constants K^.^.^ 7^ 0, contrary to what is assumed in [8]. Next, although Q 1 [y] 

' ' 71 — 1 

is complex for y > 1, it can be explicitly verified both in the near horizon region and in the large z 
limit, the explicit procedure for which is discussed in the next subsection, that the imaginary part in 
this term cancels for our solution whereas it does not cancel for the solution quoted in [8] . Furthermore, 
the authors found no order e^{r^/i) term in Eq. (j4.16p whereas we find that K^.^ 1 7^ 0. Despite these 
discrepancies in the near-horizon scalar field, the stress energy tensor is unaffected at the perturbative 
order we are probing. Similarly, the First Law is also satisfied to the appropriate perturbative order for 
both our solutions and those of Ref. [S], so these discrepancies can only be distinguished at higher orders; 
we postpone a discussion of the thermodynamics of our solutions until ^ Finally, although our solutions 
for Z**^ and 5*" at order e'^{r^/i)'^ look different than the result in [S], it can be verified that our solution 
for n = 3 agrees term by term in the near horizon expansion. 



4.2.3 Matching Region 

A crucial and often difficult step in this analysis involves matching the near region solution to the far 
region solution in order to ensure a valid solution everywhere. The heuristic procedure is as follows: take 
a small-r expansion of the outer region fields, i.e. expand F°^^ around r = 0, and match these at each 
order in e and r_|_ to the large-z expansion of the inner fields, i.e. an expansion of -F*" around z = 00. 
Such a matching takes place in an area between the two regions, where r_|_ ^ r ^ 00. This then fixes the 
remaining arbitrary constants of integration in the fields as well as fixes the expansion parameters for u. 
When taking the large z limit of the inner region fields, one must first pick an order in e, take the large 
z limit at that order, transform the limit back to the original radial coordinate, r = r^z/i, and finally 
expand the result in powers of r^/i. In principle, this is a straightforward procedure; in practice it can 
be subtle and difficult. This is best illustrated with an explicit example. 

Consider the large z expansion of the order term of g^^{z) in n = 3. The difficulty arises because 
we must take the large z limit of a function that is defined via an integral of the form 

J f{x)dx (4.19) 

for some function f{x). Now, if (j4.19p converges as we send 2 — )• 00 then this integral is easy enough to 
perturbatively compute. Indeed 

fz/i j-OO POO 

/ f{x)dx= / f{x)dx- / f{x)dx (4.20) 
Jl Jl Jz/e 

where now f{x)dx is a constant that we can easily compute, albeit numerically, and the second term 
can be evaluated by taking a Taylor series expansion of /(x) near 00. A problem arises, however, if the 
integral ()4.19p does not converge. This is true of the integral appearing in our example and is, in fact, 
generically true for the majority of the inner region fields above defined through integrals. We therefore 
must find a way to correctly take the large z limit of these integrals. 

Motivated by the above discussion for convergent integrals, if f{x)dx diverges, we define Div{x) as 
the sum of the terms which cause the integral to diverge. That is, we Taylor expand f{x) at infinity and 
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collect the terms whose integral diverges and call that Div{x). Then, by definition, {f{x) — Div{x))dx 
converges and the above consideration is valid. This allows us to rewrite our original integral as 



f{x)dx 



(^f{x) — Div{x))dx — / (^f{x) — Div{x))dx + / Div{x)dx 



z/£ 



(4.21) 



Now {f{x) — Div{x))dx is a numerical constant, the second integral can be evaluated by Taylor 
expanding f{x) — Div{x) near infinity and since Div{x) is a simple power series in x, the last integral has 
an exact analytic expression. Note that we are free to include extra convergent terms in the definition of 
Div{x) since such terms will consequently disappear from to the first two integrals of ()4.2ip . not changing 
the end result. 

Applying this technique to the order term of g'^^{z), we identify f{x) = (^{Pi[1x^ — 1])'^ 
Expanding f[x) around infinity, we find 



Div{x) 



Ax 2 
3 2>x 



2x3 



(/(x) - Div{x))dx = -10.3749804991685. 



Ci 



up to 15 digits of precision. We now have the large z limit of f{x) since the other integrals are straight- 
forward to compute. Let us now detail how the matching condition is employed. 
We must take the large z limit of the entire order at e^. That is, we must expand 




f{x)dx + Kl.,+g^^^^ 



(4.22) 



at large z. Next, we transform back to the coordinate r = zr^jt and series expand in powers of Tj^ji. 
The result for ()4.22p up to order {r+/tf is 



in 
5.2 



2r' 



3£2 £2 



19 



24 



3:2,2 



2, 

3 log 



2^ 
+ 3 log 



(4.23) 



This is now to be matched with the small r limit of the order term of g°^^{r). Direct calculation using 
the fields in Appendix B (with the constant of integration, C2, not yet fixed) yields 



out 
5e2 



2 2r2 rl 

■3 + ^+ £2" 



' 4r4 



+ + Y + 3'^2 log 



83 2 , 

1 — log 

9 3 ^ 



+ 0(r' 



(4.24) 

The order {r^/£)^ terms cancel so we continue to the order (r+/£)2 terms: Eq. (|4.24p has a term of order 
while Eq. (j4.23p does not, meaning that we must set C2 = 0. The log[r/£] terms cancel and we are 
left with a condition on the constant ^^3.2 2- Explicit calculation reveals 



-'^3;2,2 



607 ^ 2 ^ 

-^^ + 3^°^ 



(4.25) 



which has a numerical value of 1.94442... + | log • This result differs from the value of |g| + 1 log[2] + 
I log [^] ~ 1.86642... + I log [^] in Ref. [8] well beyond the margins of numerical error. This discrep- 
ancy arises because in Ref. [8] an approximation scheme was used to evaluate the large-z limit of these 
integrals [T6|, whereas our results are exact. 

Following the procedure outlined above, in n = 5, 7, 9 we find the corresponding constants to be 



^5-2 2 = 9.020251036253986, 



K|.2 2 = 8.747040877884002, 
Similar considerations show that K^.-^^ 3 are given by 



-"■9;2,2 



9.031578102169204 (4.26) 



-"^3:1,: 



14.36212300918522 



2l7r , 
^^log 



-"■5;!, 3 



-5.48986673419516 



(4.27) 
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- 7.87633500695327 



-9.11278197442749 



(4.28) 



Evaluation of the constants Kn-2 2 requires further subtle considerations since the corresponding field 
is defined via a series of nested integrals. The procedure outlined above fails under such conditions. This 
is remedied by treating the inner integrals first by Taylor expanding the integrand at infinity and near 
the horizon, while keeping a sufficient number of terms in each series such that the two match-up at 
some intermediate radius. This integral can then be evaluated directly, with the answer then becoming 
the integrand of the next integral. This integrand can then be expanded in the same way and its answer 
used as the integrand of the next integral in the nest. Finally, the last integral in the nest can be treated 
as discussed above to obtain the large z limit of the field. This rather tedious procedure yields 



K^.^ 2 = 64.6979... - 10 log ^ 2 = 306.637... (4.29) 



5;2,2 

K^2,2 = 958.528... K^^2,2 = 2451.23.. 



where we have kept only 6 digits of precision due to the numerical error introduced by the repeated 
matching of Taylor series. We again have a discrepancy well beyond the bounds of numerical error 
with the results of [8], which found Kf^2,2 = 87.5209... - 10 log [^] , again because of the approximation 
procedure used [16]. 



5 Thermodynamics and Physical Properties 

Since our solution is invariant under the Killing vector field (|2.7p . it must satisfy the First Law of 
Thermodynamics, which follows from a Hamiltonian derivation of the first law [T7j- In this section we 
work out the thermodynamic quantities of our hairy black hole solution and verify that the first law holds 
in each dimension up to the appropriate perturbative order. 



5.1 Asymptotic Charges 

In an asymptotically anti-de Sitter spacetime {M,gab) there is an ambiguity in defining "the" asymptotic 
metric because certain metric functions diverge and thus there does not exist a smooth limit to the 
boundary. To get around this subtlety, Penrose proposed a conformal completion {M,gab), where gab = 
Cl'^gab for some conformal factor Cl, such that the boundary of M is reached by the smooth limit J7 — )• 0. 
In particular, by virtue of Q vanishing (smoothly) on the boundary, Vil can be used as a radial direction 
near infinity and the subtleties of taking infinite distance limits in the physical metric gab reduces to local 
differential geometry of fields in the conformal completion gab over finite distances. Using the conformal 
completion with reflective boundary conditions, Ashtekar and Das elucidated how to properly deflne 
conserved charges in asymptotically anti-de Sitter spacetimes [18] , a procedure that was extended to the 
rotating case in |19j . Furthermore, it has been pointed out that while there exist results in the literature 
in disagreement with the Ashtekar-Das formalism, such definitions of mass and angular momenta fail to 
satisfy the first law |20j . We therefore restrict our attention to the Ashtekar-Das formalism for computing 
the mass and angular momentum, briefiy detailing the procedure. 

We start by introducing the conformal metric gab = ^^Qab where the physical metric is given in ()2.2p 
and the appropriate conformal factor is Vt = 1/r. The conformal metric explicitly looks like 

= -fgdt^ + + h{dx + Aidx^ - ^Idtf + gijdx'dx^ (5.1) 

where g, h and Q are our previous (physical) metric functions written in terms of the conformal radial 
coordinate, O, and / = O^/ now admits a smooth, finite limit to the conformal boundary, X, defined by 
= 0. The Weyl tensor, Cabcd-, of the conformal metric (|5.ip vanishes on X but ^abcd = lim ^^"^^afecd 
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admits a smooth limit, the electric part of which, £ab, is used to define conserved quantities. Using the 
transformation properties of the Weyl and Ricci tensors under conformal rescalings and imposing the 
Einstein equations, it can be shown that Sab satisfies the continuity equation [181 HH] 



D^Smd = -8n{D - 3)Tabn''h''m (5.2) 

where Ha = is the normal vector to I, Da is the derivative operator compatible with the induced 
metric hab on X and T^;, = l\mtl^~^Tab admits a smooth limit where Tab is the matter stress tensor 

in Eq. (|2.6|) . If there is no net flux of the matter stress tensor on I, then the corresponding charge is 
conserved. 

To set this up explicitly, the electric part of the Weyl tensor takes the form £ah=^'^^ambn'n^T^"' where 
the symbol = is used to denote equality on X. We also have that Ua is spacelike, which means X is 
timelike. Consider now a t=constant slice of X: this is a -D — 2 dimensional spacelike hypersurface, S, 
with induced metric 7 and normal vector ta = dat. For any conformal Killing vector, the continuity 
equation (|5.2p leads to the conserved charge 

= ± [ Sabet'dE (5.3) 

87r(n - 1) 

where the it sign corresponds to a timelike/spacelike conformal Killing vector respectively and n = D — 2 
is not to be confused with n^. 

Although our solution is invariant under a single Killing field, both dt and are asymptotic Killing 
vectors since the scalar field vanishes asymptotically. Note that the vanishing of the scalar field at 
infinity also ensures there will be no net flux of matter fields on the conformal boundary, and so the 
charges defined above are strictly conserved, dt and are then the desired conformal Killing vectors 
and they lead to a conserved energy and angular momentum respectively. A direct calculation of the 
conserved charges yields 

(n + l)7r'^r-^ 
16 (^)! 



^" 16 (^)! 

where Cf,Ch and Cq are the leading order corrections to the asymptotic fields /, h and appearing in 
the boundary conditions ()3.3p . Using the far region field expansions ()4.6p . ()4.8p and ()4.9p . along with 
the catalogued fields in Appendix B, these coefficients are easily obtained. The resulting perturbative 
expansions of the energy and angular momentum charges for n = 3, 5, 7, 9 are catalogued in Appendix A. 



5.2 Near Horizon Quantities 

Boson stars are horizonless geometries so their thermodynamics are completely governed by the above 
subsection. For the case of black holes, the presence of the horizon introduces a temperature and an 
entropy which also enter the first law. We now wish to find these quantities. 

The norm of the Killing vector ()2.7p is = —fg + r^(a; — fi)^, which is null at the horizon by 
virtue of the condition Qh = ^ and the fact that / vanishes at the horizon. The Killing vector under 
which our solution is invariant is therefore tangent to the generators of the horizon. The event horizon 
is therefore also a Killing horizon and thus has a temperature Th = For a metric of the form (j2.2p . 
a straightforward calculation yields 

Th = -^/'(r+)7^ (5.6) 
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Furthermore, any solution to Einstein's equations with or without cosmological constant obeys the 
Bekenstein-Hawking area-entropy law. Spatial sections of the horizon have an induced metric of a 
squashed n-sphere 

dsff = r^{h{dx + Aidx^f + gijdx'dx^) (5.7) 

so the entropy takes the form 



S=^rl^JW+) (5.8) 

where is the area of a (round) unit n-sphere. Now, the near region field expansions ()4.12p . ()4.14p and 
()4.15p . along with the fields catalogued in Appendix B yield the perturbative expansions of (j5.6p and 
(|5.8p . These expansions for the entropy are also catalogued in Appendix A for n = 3, 5, 7, 9. 

With the thermodynamic charges and potentials catalogued in Appendix A, it is straightforward to 
verify that the first law holds up to the appropriate perturbative order in each dimension. For boson 
stars, obtained in the r+ — t- limit, the first law takes the form 

d^ = wdJ (5.9) 

while for the hairy black holes, the first law takes the form 

dE = THdS + udJ. (5.10) 

Upon examining the thermodynamic quantities in Appendix A, we see that the angular momentum 
is primarily carried by the scalar field. Furthermore, the scalar field carries more and more of the angular 
momentum as the space-time dimension increases: the leading terms in the perturbative expansions are 
due to the scalar field at orders £^,...,6"""^ and the next terms are a mixture of scalar field and black 
hole at orders r^e'^ where p + q = n + \. Similarly, for n = 3 the energy has contributions at the same 
perturbative order from the black hole and the scalar field since the leading terms in E^, are of order 
and r\. However, for n > 3 the majority of the energy is carried again by the scalar field: the leading 
terms in En are at orders e^,...,e"'~^ and the next terms are at orders ^"""^ and e"~^. This is a rather 
striking feature of the perturbative regime of these hairy black holes: in n = 3 the thermodynamics are 
governed both by the scalar field and the black hole while in n > 3 the thermodynamics are dominated 
by the scalar field. 



6 Discussion 

We have constructed perturbative solutions describing asymptotically AdS rotating hairy black holes and 
boson stars in dimensions D = 5,7,9,11. Apart from several technical discrepancies, we are in general 
agreement with previous results [8] for the D = 5 case. Such solutions describe lumpy massless scalar 
hair co-rotating with a black hole and are invariant under a single Killing field. This is made possible by 
a particular choice of scalar field ansatz, whose stress tensor shares the same symmetries as the metric. 
These are the first known examples of asymptotically AdS black holes that are invariant under one Killing 
vector; all previous AdS black holes were stationary and hence had at least two Killing vector fields. 

The hairy black hole solutions constructed herein describe a two-parameter family of solutions, char- 
acterised by r+ and e, which is connected to the two-parameter family of Myers- Perry- AdS black holes, 
characterised by r_|_ and a; the introduction of the scalar field and requiring invariance under a single 
Killing vector fixes the hairy black hole angular momentum, leading to a two-parameter family instead 
of three. The Myers-Perry- AdS family is obtained in the limit of vanishing scalar hair. As pointed out in 
[S], sufficiently low frequency perturbations in the Myers-Perry- AdS solution will lead to a super-radiant 
instability where the amplitude of the perturbation will grow at the expense of the black hole rotation 
energy. During the instability, there are multiple frequency parameters corresponding to the various 
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perturbation modes and the solution therefore is not invariant under any KiUing fields. The end result 
of this instability, at least in the perturbative regime, will only have the lowest frequency mode present 
since higher frequency modes will get eaten by the horizon as the lower frequency modes continue to mine 
rotational energy via super-radiance. The hairy black hole solution considered herein, then, is the end 
result of a super-radiant instability of the odd-dimensional Myers-Perry-AdS black hole; the condition 
that the angular velocity of the scalar field matches the angular velocity of the horizon is a statement that 
the mode has extracted all the rotational energy it can and super-radiant scattering is no longer possible 
for that mode. Furthermore, these considerations motivate our only considering the ground state radial 
mode (|3.8p of the scalar hair and ignoring all higher radial modes since this corresponds to the lowest 
frequency mode. 

Our two-parameter family of solutions further admits a one-parameter family limit that describes 
rotating boson stars, which are horizonless geometries with a rotating clump of scalar field condensate. 
This (continuous) limit corresponds to r+ — )■ 0, which means that a perturbative hairy black hole is 
interpreted as a small black hole being added to the center of a boson star. This interpretation is 
explicitly clear by the double field expansion in powers of e and r+. One can first perform an expansion 
in powers of e to obtain the boson star solution, then expand in powers of r+ to obtain the hairy black 
hole: the second expansion corresponds to adding a black hole to the center of the boson star. It is 
possible to add a black hole to the center of a rotating boson star only if w = Oi^, which is explicitly 
demanded by the equations of motion, so that there is no net flux of scalar field across the horizon. 

The present work is a step in the direction of further understanding AdS solutions with one Killing 
vector, but it is far from complete. We have only considered the perturbative regime where energy 
and angular momenta are small by construction. In [H], non-perturbative results were also investigated 
in D = 5 and a rather rich and interesting thermodynamic structure was unveiled. The drawback is 
that one must construct solutions numerically because of the highly non-linear and coupled equations of 
motion. The benefit is that the physics is much more interesting; we refer the reader to Ref. [8j for a full 
discussion. As noted in ^ there are terms in the equations of motion which are absent in D = 5 and it 
is unclear from our perturbative analysis whether these terms have any interesting physical significance 
to higher dimensional hairy black holes. The only way to uncover this potential new physics would be to 
numerically construct hairy black holes in D = 7,9, 11 and compare them to the numerical solutions in 
D = 5. This is not a light undertaking and it is currently left for future considerations. 

It is also unclear from the present construction whether it is possible to construct asymptotically 
AdS hairy black holes with one Killing field in three-dimensions; such a solution would correspond to a 
BTZ black hole with lumpy scalar hair. This would be a desirable solution to have since L> = 3 is the 
quintessential playground for investigations of quantum gravity. However, since there does not exist a 
non-trivial Hopf-fibration of the one-sphere the prescription used in this paper is inappropriate to search 
for solutions in three-dimensions. This is related to the D > 5 Myers-Perry black holes being disjoint 
from the BTZ black holes. One must then use a different technique to construct three-dimensional hairy 
black holes, which are presumably analogously connected to the BTZ black hole in the limit of vanishing 
scalar hair. Likewise, there is nothing physically preventing analogous even-dimensional hairy black hole 
solutions from existing so it should be possible to repeat this analysis in even dimensions. The hurdle in 
this respect will be finding the appropriate form of the scalar fields such that the stress tensor shares the 
symmetries of the metric. 

Finally, it would be desirable to study these hairy black hole solutions from an AdS/CFT perspective 
to see what they correspond to in the dual gauge theory. The interpretation of our solutions from a 
boundary gauge theory perspective is currently an open issue and certainly warrants future investigation. 
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A Conserved Charges and Thermodynamic Quantities 



In this Appendix we catalogue the thermodynamic charges and potentials entering the first law in space- 
time dimension D = n + 2, for n = 3, 5, 7, 9. The boson star quantities are obtained by the limit r+ — )• 0, 
except for the temperature, which exhibits the usual Schwarzschild-like divergence as r+ — t- 0; the boson 
star temperature is zero. The constants 2 appearing in the expressions for the temperatures are 
calculated in ^4.2.3[ 
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B Perturbative Fields 

In this Appendix we catalogue all of the gravitational and scalar fields for the perturbative boson stars 
and hairy black holes in spacetime dimension D = n + 2 for n = 3, 5, 7, 9. The fields are labeled as Fn-p^g, 
where p denotes the order in e and q denotes the order in r_|_ . Note that Li2 [x] is the dilogarithm function 
while Pu[x] is the Legendre function of the first kind. 
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+ 210(3428325043311168007r2 + 1261851310575329737)r22£i6 _^ 42(23998275303178176007r2 
+ 4862510402105676583)r2°^^^ + 4158(242406821244224007r2 + 13837058565923907)r^® 
+ 10(71994825909534528007r2 - 6011370596731037917)ri6£22 _^ io(35997412954767264007r 

- 8283514692990005221)r^'^^2^ + 190(631533560609952007r2 - 250880065213913363)r^2£2^ 
+ 128(187486525806079507r2 - 116667459449721109)ri°^28 _^ 7ii36(30668879206007r2 

- 33694272545267)r^£^° - 163111679979258465r'^£32 - 18207747905205495r^£^'^ 

„2n r£2(^2^^2) 



'£2 



-log 



1+^ 



log 



+ 4102358400 ^^ 



1 + ^ 



I 10014327904r^° 



+ 225322377840r2^£2 + 2453510336480r2''^'^ + 17349823093680r2'^£^ + 91480885403040r2^^* 
+ 360803126340498r20£io + 1040279830382820ri^^^2 _^ 2178047537719290r^6£i'^ 
+ 3321710525234040r^^^^^ + 3698150956946320r^2£^*^ + 2991052821746808r^°^2° 
+ 1729094535021140r^£^2 + 690952595079480r^£^'^ + 179411207692455r^^2^ 

+ 26968400204190r2^2'^ + 1682848959147£^° 



/9.40 = -. ^ f 937254465799230r32 + 19681794632595465r3°£2 + 196811905684882635r2®£^ 

859498662576153600 V 

+ 1246433118183419250r2''^^ + 5608737609387866100r2^£^ + 19068904767424154265r22£^° 
+ 50684524900929500343r2°£^2 + 107062480455455025288r^*£^^ + 180579919814518993740r^^£^^ 
+ 242609187466772524790r^'^£^^ + 257756682509308369842r^2^2° + 213952704807502670100r^°^22 
+ 135934002405968654888r*£24 _^ 63513374918908684770r^^2^ + 19770363913882326030r^^2^ 

+ 2720901281792621928r2£30 + 133340623834666248£32 

09.40 = ^ ( 549149188365r2^ + 12081282144030r^2£2 + 126853462512315r3°£'' 

859498662576153600 V 

+ 845689750082100r2^£^ + 4017026312889975r2''^^ + 14460723267227622r2^^^° 
+ 219306855955192929r22£^2 + 2041707149885182032r2°^^^ + 986674732 1987487030r^^£^'' 
+ 29286912971461691700r^6£^^ + 58562593998521246634r^'^£2° + 82904154261107854584r^^^22 
+ 85422451772928300030r^°^2^ + 65138644777217991980r*£2^ + 37242538451355126890r^£^^ 

+ 11936492827687331904r^^^° + 1768012312961670776r2£32 + 84495843563523392^^^ 

/19.40 = ^- 1 61016576485r30 + 1281348106185r28£2 + I2813481061850r26£4 

' ' 95499851397350400 V 
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+ 81152046725050r2'^£^ + 365184210262725r22£^ + 1241596238094513r2°^^° 
+ 3310396957607208r^*^^2 _^ 7089196246488360r^6£i'^ + 12381092092392030r^^^^^ 
+ 17783794333160310r^2£^^ + 21034536810890148r^°£2° + 20300574058072308r^£2^ 
+ 15541844723695290r^£2^ + 8797952237444370r^£2^ + 3015895473831240r2^2^ 

+ 14972530721 1480£^° 

1^9-4 = ( 5440474418585r^^ + 114249962790285r^°^2 + 1142499627902850r2^£^ 

' ' 59687407123344000 V 

+ 7235830976718050r2^^^ + 32561239395231225r^'^£* + 110708213943786165r22^^° 
+ 295268817053306040r^°£^2 + 632417867298394200r^*^£^'^ + 1101647489722639650r^''£^^ 
+ 1565834578187766250r^^^^^ + 1807245829726785084ri2£20 _^ i674i73405090497584ri°£22 

+ 1221884704722504550r^£^'^ + 681865813132481250r^^2'' + 273777806813998500r'^£^^ 
+ 64940197030107672r2^^° + 6598585567818972£32 

Hg-so = I 196278755076119819922240261840r^2 

12928610477426464084962297600000 V 

+ 4302621620873952400266415533360r^°£2 + 45035990011802132234460421933080r^^£^ 
+ 299476749663332693852942194632980r^'^^^ + 1419986713892478010760754706555122r^^£^ 
+ 5107957999464578930096857586903749r32£i0 _^ i4480110437447756571401543030317327r^°£^2 
+ 33171546764185616749586054262623220r28£i4 _^ 62477592200815789501550742105430020r26£i6 

+ 97921202768932011707091446069275799r^^£^^ + 128775176906725704442407617007965163r22£2° 

+ 142873004760249240821838450887004074r2°^22 + 134129087251912145628690383503063720r^^^2^ 
+ 106606113471306396324075091246911420r^6^26 ^ 7i564603098088232189439094493589402r^^£28 

+ 40318687255770901091487447290055224r^2£^° + 18828704309439914264745324330571272r^°^32 
+ 7126896792786249923435515161591030r^£3'' + 2095459575089784573411052833767630r^£3'^ 
+ 438680016278348277072828543787740/^^^ + 54960718928326817581432660068800r2/° 

+ 2963660262582405195600812804780^"^^ 
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